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We introduce a new picture of vacuum decay which, in contrast to existing semiclassical tech- 
niques, provides a real-time description and does not rely on classically-forbidden tunneling paths. 
Using lattice simulations, we observe vacuum decay via bubble formation by generating realizations 
of vacuum fluctuations and evolving with the classical equations of motion. The decay rate ob- 
tained from an ensemble of simulations is in excellent agreement with existing techniques. Future 
applications include bubble correlation functions, fast decay rates, and decay of non-vacuum states. 


Introduction — Field theories with complex interac- 
tion potentials arise in a diverse range of physical appli- 
cations: from cosmological theories of the multiverse in- 
spired by string theory [1] [2], to cold atom Bose-Einstein 
condensates [3H5], to the dynamics of protein folding [6]. 
Such theories generically possess many local potential 
minima. At the homogeneous classical level such local 
potential minima are stable. However, the uncertainty 
principle ensures that spatially-inhomogeneous quantum 
fluctuations are present, resulting in the metastablility 
of false vacuum field configurations. The decay of false 
vacuum states is expected to proceed analogously to first- 
order phase transitions driven by statistical fluctuations 
in condensed matter systems [7][8], with bubbles of a new 
phase nucleating then expanding into the ambient false 
vacuum. The prevailing, if implicit, wisdom is that bub- 
bles nucleate via non-classical paths in field space, analo- 
gous to tunneling through a barrier in quantum mechan- 
ics [9HI6]. We propose an alternative picture of vacuum 
decay, in which the classical evolution of the field from 
some initial realization of the false vacuum fluctuations 
leads to the emergence of bubbles. As convenient nomen- 
clature, we will refer to these as classically-forbidden and 
classically-allowed nucleations, respectively. 


We re-examine the nature of bubble nucleations dur- 
ing false vacuum decay in relativistic scalar field theory 
using semiclassical scalar field lattice simulations. In our 
framework, quantum effects are included by sampling re- 
alizations of the initial quantum state of the field, then 
time-evolving using the classical equations of motion. As 
we demonstrate in the supplementary material, this cap- 
tures the time-evolution of the wave functional to O(h) 
(see also Ref. [I7] for a similar derivation). This treat- 
ment is analogous to the truncated Wigner approxima- 
tion [18] often employed in atomic physics. We ex- 
plicitly demonstrate the existence of classically-allowed 
paths connecting these initial realizations of the quan- 
tum false vacuum, to a subsequent state with nucleated 


bubbles of “true vacuum”. This demonstrates that false 
vacuum decay can proceed via classically-allowed dynam- 
ical evolution. Furthermore, we compare the nucleation 
rates inferred using classical statistics for the initial fluc- 
tuations (i.e.,ignoring quantum interference), and find 
good agreement with instanton-based predictions which 
cannot tackle the time-dependent nature of the prob- 
lem. This suggests that the decay of the false vacuum 
proceeds via classically-allowed bubble nucleation events. 
The supplementary material shows how this approach 
can be interpreted as an expansion of the full quantum 
dynamics in h. 

In this work, we focus on the particular field theory 
emerging from the cold atom analog false vacuum pro- 
posal of Fialko et al [B}H5). However, we do not expect 
the presence of classically-allowed decay channels to de- 
pend strongly on the specific choice of field theory. The 
possible connection to experiment is particularly exciting 
in this context, as it provides a window to experimentally 
confirm the dynamics of vacuum decay, and possibly un- 
cover additional effects arising from quantum interference 
such as Anderson localization [20} [21]. 

Semi-Classical Statistical Formalism — We study 
relativistic false vacuum decay using semi-classical sta- 
tistical lattice simulations. Throughout we work in units 
with A = c = 1. In order to generate sufficient statis- 
tics to compute decay rates, we will work with fields in 
1 + 1-dimensions, although we do not expect considera- 
tion of higher spatial dimensions to qualitatively change 
our conclusions. 

The fields are initialized as 


(x,t = 0) = dy + 64(x) 


where the fluctuations 5¢ and 6¢@ are realizations of 
random fields with the same statistics as the quantum 
Minkowski vacuum fluctuations, and @, is the mean field 
value in the false vacuum. Using these as initial condi- 


tions, we solve Hamilton’s equations 


dp, do 2 ; 

Ho? and a=-Vvb-Vig) (2) 
with a tenth-order accurate Gauss-Legendre time inte- 
erator [23]. The spatial Laplacian is computed using 
a Fourier collocation stencil, resulting in periodic bound- 
ary conditions. This combination conserves energy to 
O(10~'°) with a maximal pointwise error of O(10~'%) 
as either the grid spacing or time-step is varied. Addi- 
tionally, our time-integrator is symplectic so that phase 
space volume is preserved. Finally, we verified that our 
numerical simulations are time-reversible. First we in- 
tegrate forward in time well into the regime where the 
false vacuum has decayed completely. By reversing the 
flow of time and backward integrating for an equal time 
interval, we are able to recover the initial state of @ and 
@ to a precision of O(10~'>). Further, if we continue to 
backward integrate in time, we again see the initial state 
decay since the false vacuum is not a true eigenstate of 
the Hamiltonian. 

This approach does not capture quantum interference 
or the contributions of non-classical paths, which as 
we demonstrate in the supplementary material enter at 
O(h?) in our approach. Nonetheless, there are many 
cases in which this semiclassical approach gives an ac- 
curate description of the statistics of the field(s). Some 
well-known cases include post-inflation preheating dy- 
namics and relativistic heavy ion collisions [24}{28]. In 
these standard scenarios, individual Fourier modes be- 
come highly squeezed while in the linear regime, thus 
entering the classical wave limit and justifying a classical 
treatment of the field dynamics and statistics. The case 
considered here is novel, since once a bubble has formed 
in the medium the fluctuations are no longer statistically 
homogeneous, and the growth of modes leading to its for- 
mation may not occur in the linear regime. However, a 
bubble is a highly classical field configuration, so we ex- 
pect that many aspects will be properly captured by our 
simulations. Regardless of how accurately these semiclas- 
sical simulations capture the full quantum statistics, they 
are more than sufficient to demonstrate the existence of 
classically-allowed bubble nucleation trajectories. 

Classically Allowed Decays— We now use the 
framework outlined above to study the dynamics of a 
scalar field ¢ initially trapped in a false vacuum. Mo- 
tivated by the analog false vacuum setup in cold atom 
Bose-Einstein condensate (BEC) systems [3H5], we con- 
sider the following action 


P _ (Vey ov _2¢ 
5 5 va ( cos a 5 sin? a (3) 


In the BEC experiments, \ is an experimentally tunable 
parameter controlling the depth of the false-vacuum well. 
For \ > 1, this theory posesses an infinite sequence of 


i= 





false vacuum extrema at ¢ = (2k + 1)n with k € Z. 
Instead, for A < 1 the false vacua are local maxima of 
the potential, and states centered around these points 
undergo spinodal decay. 

Here, we approximate the initial vacuum fluctuations 
as those of a free massive scalar field with mass squared 
m2, = V"(¢ = 1¢0) = Vobg?(-1 + 2). We sample 
fluctuations of the field 6¢ and its time derivative 5¢ as 
uncorrelated Gaussian random fields with spectra 


(Oidu) = Ahk) (Bide) = SHak—R), (4) 


where wz = k? + m2— and 6(k — k’) is the Dirac delta 
function. 

It is convenient to work with the dimensionless variable 
@ = $/¢0, along with dimensionless units £ = pt and 
Z = px where yp is a constant with dimensions of mass. 
Explicitly, for a one-dimensional discrete lattice of finite 
size L, 


“ 11S, a 
5¢(x) = ——= ) ee *__" __, 65 
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where Gy are realizations of complex Gaussian random 
deviates with unit variance, which we generate using 
Qn = J—Inane!2™= where G, and 6, are uniform ran- 
dom deviates on the unit interval. We have included a 
sharp cutoff in the spectrum, which must be less than or 
equal to the Nyquist frequency and the wavevectors are 
given by ky = 20 The 6¢ fluctuations are initialized 
analogously. 

Figure [I] shows an example classical field evolution in 
one spatial dimension for a single initial field realization 
described above. We see a bubble (i.e.,domain wall- 
antiwall pair) nucleate from the initial vacuum fluctu- 
ations, then subsequently expand into the surrounding 
false vacuum. To our knowledge, this is the first explicit 
demonstration that bubble nucleation in false vacuum 
decay can occur via classically-allowed field evolutions in 
pure relativistic field theory. Aside from the nature of 
quantum vs. classical statistics, this is analogous to bub- 
ble formation in thermal [29] (or other classical statis- 
tical [30]) phase transitions. Thermal nucleation events 
are also typically modeled as random, but at a micro- 
scopic level the full set of field fluctuations (including the 
“thermal bath”) are simply following the classical equa- 
tions of motion. Bubble nucleations are special classical 
trajectories through phase space that describe the coher- 
ent creation of a bubble. In the quantum case we are 
modeling, the fluctuations instead descend directly from 
the initial state of the field rather than an external bath. 


Let us now consider how to interpret the results of 
our simulations. A complete time-dependent description 





FIG. 1. A classically-allowed false vacuum decay for the La- 


grangian density with A = 1.2 and ¢o = ree The initial 


fluctuations 6¢ and 50 are realizations of Gaussian random 
fields with spectra mimicking Minkowski false vacuum fluc- 
tuations. We observe the dynamical emergence of a bubble 
with field localized near a true vacuum (cos ¢ = 1) embedded 
within the false vacuum (cos ¢ = —1). This explicitly demon- 
strates the existence of classically-allowed paths for false vac- 
uum decay. 


requires obtaining the wave functional 


(9(x)|¥()) = [ D0.(olO (lt) o)(6s(20))¥Ct0) (6) 


for the full scalar field configuration as a function of time, 
given the initial wave-functional |Y(to)). However, any 
observer watching the false vacuum boil will not have ac- 
cess to the full quantum state of the field as encoded in 
|W(t)). Rather than seeing a quantum superposition of 
nucleation events, such an observer will see N, bubbles 
nucleate with spacetime locations {(x;,t;),i = 1, Np}. A 
priori, we can only predict probabilities for N, and these 
spacetime positions. Given that an actual set of nucle- 
ation events has taken place, the observer should be able 
to describe the field evolution subject to the constraint 
of these observed nucleations. Since the observer will 
not have perfect knowledge of the multi-bubble state, 
their observations should be marginalized over fluctua- 
tions around the constrained field path. The result pre- 
sented here shows that there exist constrained paths de- 
scribing purely classical time-evolution connecting final 
multi-bubble states directly to realizations of the initial 
vacuum fluctuations. Since these classical paths extrem- 
ize the action, we expect that they will provide an im- 
portant contribution to the propagator determining the 
evolution of the wave-functional. 

To contrast our approach, let us instead briefly con- 
sider the standard Euclidean (i.e.,instanton) formal- 
ism [9] [10| [14]. We Wick-rotate to Euclidean time, then 
search for extremal paths of the Euclidean action con- 
necting the false vacuum at infinity to the true vac- 
uum at the origin. These extremal paths are known 


as bounce solutions. Nearly all existing treatments as- 
sume the bounce is spherically-symmetric in Euclidean 
signature. This dramatically reduces the phase space of 
allowed field configurations, simplifying it to a quantum 
mechanics problem. Nearly all existing treatments as- 
sume the bounce is spherically-symmetric. In Lorentzian 
signature, the bounce describes a spherically-symmetric 
bubble that contracts from infinite size in the distant 
past, reaches a minimum turnaround radius, and then 
re-expands. The nucleation event itself is modeled by 
artificially removing the contracting phase, then pasting 
the remaining expanding phase into the ambient space- 
time. The pasting step is interpreted as a classically- 
forbidden quantum tunneling process, analogous to the 
decay of an unstable particle in quantum mechanics. 

However, the tunneling interpretation is fueled by anal- 
ogy with quantum mechanics, which as mentioned above 
requires a drastic reduction in the allowed phase space 
for the field evolution. We have demonstrated that the 
neglected fluctuations can play a pivotal role around 
the time of nucleation, and allow for bubbles to nucle- 
ate via classically-allowed trajectories. Specific fluctua- 
tion configurations that may lead to bubble formation 
are: (1) local peaks or valleys in the initial fluctuations 
around the false vacuum that probe the nonlinear struc- 
ture of the potential and act as bubble nucleation seeds, 
or (2) the repeated mildly nonlinear interactions of prop- 
agating waves that eventually lead to the development of 
a local field peak from which a bubble may form. The 
latter effect means that energy within a small region may 
not be conserved due to energy transport via currents. 

False Vacuum Decay Rates— Thus far, we showed 
the existence of classically-allowed paths through field 
space describing false vacuum decay via bubble nucle- 
ation. Given such paths exist, we expect they play 
an important part in the false vacuum decay process. 
They may either be a competing process with classically- 
forbidden bubble nucleations; or, given the drastic di- 
mensional reduction used in the Euclidean formalism, 
may actually be a more complete description of the de- 
cays captured by instanton techniques. To investigate 
this issue, we estimate the bubble nucleation rate using 
our real-time simulations and compare it with the instan- 
ton predictions. 

In D = d+ 1 spacetime dimensions, the Euclidean 
prediction for the nucleation rate per unit volume is 


Te (2 - 5 (seeeselenl) 


Vi \On detd2. Sp [br] 
B\P/ 
=(F) 80000), (7) 


where B = Spidp] — Seley] is the difference in the Eu- 
clidean action evaluated on the bounce solution and in 
the false vacuum, det’ indicates the determinant with 
zero modes removed, and 67Sp is the second variation of 


the Euclidean action. The bounce @p satisfies 
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with boundary conditions ¢g(rz = oo) = dg and 
Or, O(rE = 0) = 0. In addition, ¢g(rz = 0) should 
be near the true vacuum, and the solution should pos- 
sess a single negative eigenmode to give an imaginary 
contribution to the Euclidean effective action [BI]. A 
precise calculation of the determinant prefactor is be- 
yond the scope of this paper, so we parametrize it as 
D(A, ¢0) = 9(A, ¢0) (Vod52)”, where g is a dimension- 
less function of order unity which includes counterterms 
to renormalize the determinant] For the purpose of ob- 
taining scaling relations, it is convenient to re-express 
as 
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= 24( 260) } sbcprjerrboo), (9) 


we 2 

where C(X) = f redre (Cmee” + [V(d) -— V(ew)]) 
and we have specialized to one spatial dimension. For 
our particular model, V(¢) = Vo (- cos @ + a sin? 0), 
The factor of 2 accounts for the independent instanton 
decay paths in this model. 

As is clear from Fig. |1} the emergence of a bubble is 
a “fuzzy” process which makes a precise determination 
of the nucleation time somewhat ambiguous. We there- 
fore adopt the following operational definition. Imagine 
an ensemble of observers each watching a finite region of 
the false vacuum boil. We ask each of these observers to 
provide us with a history of (cos ¢)y within their region, 
where (-)y represents a spatial average. To obtain statis- 
tics, we make one additional assumption — each field 
trajectory (and thus each observer) is weighted by the 
probability of the initial realization of the fluctuations. 
In the particular case considered here, where we assume 
the initial state is Gaussian and draw samples from an 
ensemble of Gaussian random fields, each simulation is 
weighted equally. We select a threshold, and declare that 
the observer’s region has decayed at the first passage of 
(cos ¢)y past this threshold. To reduce the effects of mul- 
tiple bubbles and fluctuations, we choose the threshold 
to be ép-+n,Acr for some constant n,. Here Gp and Acrp 
are, respectively, the ensemble average and standard de- 
viation of (cos ¢)v on the initial time slice. We then run 
ensembles of 5000 lattice simulations with varying do and 
side length LD for X = 1.2, and empirically determine the 
fraction of undecayed trajectories Fyurvive(t). We verified 
Feurvive iS insensitive to the choices of both the lattice 


1 When the bounce solution possesses a hierarchy of scales, such 
as in the thin-wall limit, g may not be order one. 
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FIG. 2. Decay rates of false vacuum initial states as ¢o is var- 
ied for our ensembles of lattice simulations (orange dots). For 
comparison, the blue solid line is the analytic prediction 

with g = —1+ A? = 0.44 to set the overall scale to the false 
vacuum mass, and the gray dashed line is the leading-order 
prediction (i.e., T = Ae~2"C)%0 with the constant A normal- 
ized so the curves match at ¢o = 1). We fixed \ = 1.2 and 
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spacing dx and the threshold n, for 5 < ng < 25. Af 
ter an initial transient that depends on the simulation 
volume L, Fyurvive enters a regime of nearly exponential 
decay. We fit this exponential tail Fyurvive = e7!~*) 
with parameters I and to to extract the decay rate and 
verified that i is independent of the simulation volume 
D. 

In Fig. [2] we compare the decay rates extracted from 
our simulations with those of the instanton formalism 
(obtained using the method of Ref. [82]). We also include 
the leading order prediction to compare to an exponen- 
tial falloff at large @9. To capture all relevant dynam- 
ical modes, we set the spectral cutoff key, = 12247 Yo 


25/2 40 
above the wavenumber where the Fourier amplitude of 


the bounce profile ¢g dropped below 107'°. In this 
regime, increasing keyt should simply renormalize the 
field theory. Given the unknown determinant prefactor, 
exponential sensitivity of the decay rate to the bounce 
action, and multi-bubble and fluctuation effects entering 
into our decay time extraction, our simulations match the 
instanton prediction very well. In particular, the ampli- 
tude is the correct order of magnitude, and the rate of 
exponential decay at large ¢g closely matches the Eu- 
clidean calculation. 

We also varied \ at fixed ¢2C(A) and found the time- 
evolution of the survival fraction to be nearly indepen- 
dent of the choice of ». This suggests not only that 
classically-allowed decay paths are a key process in false 
vacuum decay, but also that these classically-allowed 
paths may be the same ones being captured by the stan- 
dard instanton calculation. In order to address this latter 
issue, the results presented here motivate a deeper inves- 
tigation of the function g, renormalization effects, and a 





more sophisticated bubble extraction algorithm. 

Conclusions — We outlined a novel real-time frame- 
work to study quantum false vacuum decay, based on 
semi-classical stochastic lattice simulations. This pro- 
vides a step towards incorporating the full temporal evo- 
lution of a decaying region of the false vacuum. In partic- 
ular, we demonstrated that the false vacuum can decay 
via bubble nucleation along classically-allowed field tra- 
jectories, and the decay rates of such decays agree well 
with standard computations based on the instanton for- 
malism. This suggests that our formalism provides a de- 
scription of the “fuzzy” process of a bubble nucleating 
in false vacuum. We will explore how our interpretation 
generalizes to other systems of non-equilibrium bubble 
formation in a forthcoming work. 

One particularly exciting prospect is the possiblity of 
emulating the false vacuum decay process with cold atom 
systems, suggesting we may be able to experimentally 
test the picture outlined here [BH5]. In addition, our ap- 
proach opens new avenues of investigation which cannot 
be addressed within the Euclidean framework. These in- 
clude decay rates from non-vacuum initial states, inves- 
tigation of bubble-bubble correlations, and evolution in 
a regime of rapid decays. Furthermore, the emergence of 
a Classical bubble from a highly quantum initial state is 
an example of self-decoherence of a quantum field, and 
may yield insight into the process by which effectively 
classical density perturbations arise from quantum fluc- 
tuations during inflation [B3}36]. Finally, the emergence 
of a classical bubble from the quantum vacuum bears 
many similarities to the process of “measuring” a bub- 
ble nucleation event, and may shed light on the notion 
of measurement in quantum field theory and quantum 


cosmology [37]. 
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SUPPLEMENTAL MATERIAL: SEMICLASSICAL STATISTICAL SIMULATIONS AS EVOLUTION OF 
THE WIGNER FUNCTIONAL 


In this supplementary material, we demonstrate that the statistical lattice simulations employed in the main text 
is equivalent to quantum dynamics to first order in an expansion in fi. A similar derivation can be found in Ref. 
Consider the Wigner functional 


W(¢, | = [Pret setemennt (6+ “ é 








1) 
6-2) (10) 
which is normalized as 
[Perm wie. =1, (11) 


where we have absorbed some factors of 27 into the definition of DI. Using W, expectation values of an operator 
O(¢, II) are obtained from 


(06,1) = / Dept Wlo,MOw(e,1) (12) 


where Ow/(¢, II) is the Wigner symbol of O. For a fully symmetrized operator, Ow is obtained by simply replacing 
the quantum operators ¢ and II by classical functions. Thus, the Wigner functional acts as a quantum probability 
distribution on phase space, although one should bear in mind that it is not always positive definite. However, for a 
Gaussian state, the Wigner function is positive definite and the probability interpretation holds. 

Consider a theory with Hamiltonian density of the form 
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Working in the Schrodinger picture, the equation of motion for the density matrix / is 
dp lyfe. 
Bole a 
Using this, we can derive the evolution equation for the Wigner functional 
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where 04 acts only on the potential and 6/éII acts on the Wigner functional. Expanding the sine term arising from 
the potential in h, we obtain 


O a (Std dH Sd Se shih & bai 
In the above, we have used the classical equations of motion 
dé oH _ 
7H (17a) 
dil ) 
=" _724_v'td), (17b) 


dt 6 


Therefore, to O(h?) in the dynamical evolution, we can solve for W using the method of characteristics, with each 
characteristic satisfying the classical equations of motion from some initial @ and II configuration. This cleanly 
factorizes the effects of quantum mechanics into: (1) a contribution from the initial fluctuation statistics encoded 
in W[¢, 1, t = 0}, and (2) O(h?) and higher modifications to the classical dynamical evolution. Our approach thus 
corresponds to this approximation, with the initial state approximated as a Gaussian. Since the initial fluctuation 
spectra are set by commutation relations between @¢ and II, in this formalism the leading quantum correction in h 
arises from the initial form of W, rather than due to dynamical evolution. Note the crucial role the existence of 
classically-allowed false vacuum decay paths plays in the utility of this approach. 

When considering Fig. 2, note that the instanton calculation is also an expansion in h. The leading effect arises 
from the action on the instanton trajectory e~*', while the fluctuation determinant gives the first quantum correction 
(including renormalization effects). Thus, both the curves are approximations to the true tunneling rate and a careful 
calculation of the determinant prefactor could bring them into even better agreement. 


